THE CAUCHY PROBLEM FOR THE 3D NAVIER - STOKES EQUATIONS. 
NEW APPROACH TO THE SOLUTION AND ITS JUSTIFICATION. 

A. TSIONSKIY, M. TSIONSKIY * 

Abstract. Some known results regarding the Euler and Navier-Stokes equations were obtained by different authors. 
Existence and smoothness of solutions for the Navier-Stokes equations in two dimensions have been known for a long time. Leray 
[2] showed that the Navier-Stokes equations in three space dimensions have a weak solution. Scheffer [3], [4] and Shnirelman [5] 
obtained weak solution of the Euler equations with compact support in spacetime. Caffarelli, Kohn and Nirenberg [6] improved 
Scheffer's results, and F.-H. Lin [7] simplified the proof of the results of J. Leray. Many problems and conjectures about behavior 
of weak solutions of the Euler and Navier-Stokes equations are described in the books of Bertozzi and Majda [S], Constantin 
[9] or Lemari-Rieusset 1101 . 

Solutions of the Navier-Stokes and Euler equations with initial conditions (Cauchy problem) for 2D and 3D cases were 
obtained in the convergence series form by analytical iterative method using Fourier and Laplace transforms in paper I12| . 
These solutions were received in a form of infinitely differentiable functions, and that allows us to analyze all aspects of the 
problem on a much deeper level and with more details. Also such smooth solutions satisfy the conditions required in 1111 for 
the problem of Navier-Stokes equations. 

For several combinations of problem parameters numerical results were obtained and presented as graphs 1121 . 1131 . 

This paper describes detailed proof of convergence of the analitical iterative method for solution of the Cauchy problem for 
the 3D Navier - Stokes equations. The convergence is shown for wide ranges of the problem's parameters. Estimated formula 
for the border of convergence area of the iterative process in the space of system parameters is obtained. Also we have provided 
justification of the analytical iterative method solution for Cauchy problem for the 3D Navier-Stokes equations. 

1. Introduction. 

Approach to the solution of Cauchy problem for the 3D Navier-Stokes equations described in this paper is 
based on the form of differential equations in the statement of the problem and also limitations for the initial 
conditions and applied force. It grows from classic definition of function and classic methods of analysis. 

First we have moved non-linear parts of equations to the right sides and then solved Cauchy problem 
for the Navier-Stokes equations by analytical iterative method. By doing that we have obtained more and 
more precise results for the right (non-linear) parts of equations on each step of iterative process. 

For this purpose we have solved the system of linear partial differential equations with constant coef- 
ficients on every step of iterative process. We have obtained the solution using Fourier transforms for the 
space coordinates and Laplace transform for time. From theorems about application of Fourier and Laplace 
transforms for system of linear partial differential equations with constant coefficients we see that in case if 
initial conditions and applied force are smooth enough functions decreasing in infinity, then the solution of 
the system of linear partial differential equations with constant coefficients is also a smooth function. (Cor- 
responding theorems are presented in S. Bochner [14], V.P. Palamodov jT5], G.E. Shilov [16], L. Hormander 
[m, S. Mizohata [H], J.F. Treves [H]). 

In the problem statement for Navier-Stokes equations initial conditions and applied force are infinitely 
differentiable functions decreasing rapidly to zero in infinity. Hence, all these theorems are applicable in 
this case. Next, by using theorems from |16j regarding Fourier transforms for infinitely differentiable and 
decreasing to zero in infinity functions we have shown that all functions of Cauchy problem for Navier-Stokes 
equations are staying infinitely differentiable and decreasing to zero in infinity also rapidly enough in Fourier 
transforms on any step of iterative process. 

After that we have received a superior estimation for the solution of the Cauchy problem for the 3D 
Navier - Stokes equations by iterative method. The purposes of this estimation are: 

1) to show convergence of the iterative method; 

2) to obtain analytical form of the first and second steps of the iterative process; 

3) to receive estimated formula for the border of convergence region of the iterative process in the 
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space of system parameters. 

Also we have provided justification for the solution by analytical iterative method of Cauchy problem 
for the 3D Navier-Stokes equations. While doing so, we have introduced perfect spaces of functions and 
vector-functions (I.M. Gel'fand, G.E. Chilov [10]), in which we have looked for the solution of the problem. 
We have demonstrated equivalence of the solution of Cauchy problem in forms of differential and integral 
equations. As a further step we have proved existence and uniqueness of the solution of Cauchy problem in 
all time range [0,oo) by using the fixed point principle ( L.V. Kantorovich, CP. Akilov [21], V.A. Trenogin 
[22], W. Rudin [21, W.A. Kirk and B. Sims [24], A. Granas and J. Dugundji [25], J.M. Ayerbe Toledano, 
T. Dominguez Benavides, G. Lopez Acedo [55] ) . For this purpose three following theorems were proven in 
this paper: 

Theorem 1: Integral operator of the problem is a contraction operator; 

Theorem 2: Existence and uniqueness of the solution of the problem is valid for any t G [0,oo); 
Theorem 3: Solution of the problem is depending on t continuously. 

By using the concept of comparable norms we have shown that the energy of the whole process has a 
finite value for any t € [0,oo). 

2. The mathematical setup. 

The Navier-Stokes equations describe the motion of a fluid in {N = 3). We look for a viscous 
incompressible fluid fllling all of here. The Navier-Stokes equations are then given by 

(2-1) ^ + E"n|^ = '^^^k - ^ + fkix,t) {xeR^, t>0, l<k<N) 

at ^ dXn OXk 

n—1 

(2.2) divu = = {x£R^,t>0) 

^ OXn 

with initial conditions 

(2.3) u{x,0) = u°{x) {x e R^) 

Here u{x,t) = {uk{x,t)) e i?^, (1 < fc < N) — is an unknown velocity vector (A^ ~ 3), p{x,t) — 
is an unknown pressure, if{x) is a given, C°° divergence- free vector held , fk{x,t) are components of a 
given, externally applied force f{x,t), v is a. positive coefficient of the viscosity (if i/ = then (|2.ip - (|2.3p 
arc the Euler equations), and A — X]^!Li '^P' Laplacian in the space variables. Equation (12. ip is 

Newton's law for a fluid element subject. Equation (|2.2p says that the fluid is incompressible. For physically 
reasonable solutions, we accept 

(2.4) Ukix,t) , ^Oas |x|-!^oo (l<k<N, l<n<N) 

OXn 

Hence, we will restrict attention to initial conditions HP and force / that satisfy 

(2.5) I d^ii^ix) I < C„k(1+ I X \)-^ on for any a and K. 
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and 

(2.6) I d!^dff{x, t) I < C„/3A'(1+ I X I +ty^ on x [0, oo) for any a, /3 and K. 
We add (— J2n=i ^"ff^ ) both sides of the equations (|2.1I) . Then we have: 

(2.7) ^ ^ uAuk - + Mx,t)- Y^un^ {xeR'', t>0, l<k<N) 

n—l ^ 

N „ 

(2.8) divu = ^ pi = {x&R^,t>0) 

n—l 

(2.9) u{x,0) = ^/'(x) (x e i?^) 

(2.10) Wfe(a;,<)^0 , ^ ->Oas |x|->oo (l<k<N, l<n<N) 

OXn 

(2.11) I a^?/'(a;) I < C„/f (1+ I a; 1)--^ on for any a and K. 

(2.12) I d"d^f{x,t) I < Cai3K{l+ I X I +<)"-f'^ on R'^ x [0, oo) for any a,/3 and K. 

We shah solve the system of equations (|2.7p - (|2.12p by the anahtical iterative method. To do so we 
write this system of equations in the following form: 

(2.13) ^ = ,yAu,k - ^ + fA^^t) (xei?^, t>0, 1 < k < N) 

(2.14) divu, = = (xei?^,i>0) 

n—l 

(2.15) Uj(a;,0) = ?/'(a;) {x e R^) 

du 

(2.16) Ujfc(x,t)-^0 , ^Oas |x|^oo (l<k<N, l<n<N) 

(2.17) \d^ii°{x) \< CaKil+\x\)-^ on R'^ for any a and K. 
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(2.18) I dy^fix,t) I < C„/3k(1+ I X I +t)-^^ on x [0,oo) for any a,/3 and K. 

Here j is the number of the iterative process step (j — 1,2,3,...). 



N 



(2.19) fA^,t) = Mx,t) - ^^,^-_l.„^^M il<k<N) 
or the vector form 

(2.20) fj{x,t) = f{x,t) - (m,_i ■ V)uj-i 

For the first step of the iterative process (j = 1) we have: 

(uo ■ V ) Wo = 

and 

h{x,t) = f{x,t) 

3. Solution and Estimation. 

Let us assume that all operations below are valid. We will prove the validity of these operations in next 
paragraphs. 

We will solve the system of equations (|2.13p — (|2.20p by analytical iterative method with the condition 
that fjk{x,t) is known function on any step j of iterative process. Hence, we have that (|2.13l) — (|2.15p is 
the system of linear partial differential equations with constant coefficients. 

Because of that on the first stage we use Fourier transform (|A.1I) for solution of equations (|2.13p — (|2.20p 
and get: 

C/?fc(7i,72,73,0 = F[ujk{xi,X2,X3,t)] 

r d'^Ujk{xi,X2,X3,t)i f +^ r fttttii 

^[ — ^^2 J = -7sC^ife(7i>72,73,i) [use(|2.16p] 

t^°(7i,72,73) = F[ul{xi,X2,x-i)] 

Pj{ll,l2,l3,t) = F[pj {xi,X2,X3,t)] 
Fjk{ll,l2,l3,t) = F[fjk{xi,X2,X3,t)] 



k,s = 1,2,3 
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and then: 

(3-1) = _j,(^2 ^^2 _^^2^[/^.^(^^^^2^^^^^) _^ i7i^'j(7i,72,73,i) + -P,i(7i,72,73,0 

(3.2) ^^'^^^^^^J^''^^'^'* = -^(7i +72 +73)C^j2(7i,72,73,0 + «72^'i(7i, 72, 73, *) + -^^2(71, 72, 73, 

(3-3) ^^'^^^^^^J^'^^'^'' = -^^(7? +7I + 7|)t^j3(7i, 72,73,0 + «73^'i(7i,72,73,i) + -^,3 (71, 72, 73, 



(3.4) 7iC^ji(7i,72,73,0 + 72^^^2(71,72,73,0 + 73 t(?-3(7i,72,73,0 = 

(3.5) C/,i(7i,72,73,0) = [/i (71, 72, 73) 

(3.6) t^j2(7i,72,73,0) = ^72(71,72,73) 

(3.7) C/,-3(7i,72,73,0) = 1/3(71,72,73) 

Hence we have received a system of hnear ordinary differential equations with constant coefficients in 
regard to Fourier transforms p.ip — p.7p . At the same time the initial conditions are set only for Fourier 
transforms of velocity components C/ji(7i, 72, 73, 0, ^7/2(71, 72, 73, 0, C/)3(7i, 72, 73, 0- Because of that we 
can eliminate Fourier tranform for pressure Pj (71, 72, 73, t) from equations ([23J ~ p.3p on the second stage 
of solution. From here assuming that 71 ^ 0, 72 7^ 0, 73 ^ 0, we eliminate fj(7i, 72, 73, i) from equations 
(|XT|) - (Hini) and find: 

C/i2(7i,72,73,i) - — C/,i(7i,72,73,0] = -i^(7i + 72 + 73)[ ^7^2(71, 72, 73, i) - 
at 71 

72 1 r 72 1 

(3.8) C/,-i(7i,72,73,0J + [ ^"^2 (71, 72,73, i) i^ji(7i,72,73,0J 

71 71 

4[ C^j3(7i,72,73,0 - — C/ji(7i,72,73,0] = -^^(7? + 7I + 7|)[ f/j3(7i, 72, 73, - 
at 71 

(3.9) - — C/j-i(7i,72,73,0] + [ -F^j^ (71,72,73,*) - — i^ji(7i,72,73,0] 

7i 71 



(3.10) 



71^7,1(71,72,73,0 + 72^7^2(71,72,73,0 + 73 C7(-3(7i,72,73,0 = 
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(3.11) i7ii(7i,72,73,0) = C/i°(7i,72,73) 

(3.12) C/,-2(7i,72,73,0) = [/2(7i,72,73) 

(3.13) t/j-3(7i,72,73,0) = [73(71,72,73) 

On the third stage we can use Laplace transform (jA.2p , (|A.3p for a system of hnear ordinary differential 
equations with constant coefficients p.8p — (|3.10p and have in result a system of linear algebraic equations 
with constant coefficients: 

Ufi.ilul2,l3,v) = i[C7,-fc(7i,72,73,i)] k=f,2,3 



C7|(7i,72,73,?7) - — (71, 72,73, - [ f7,-2(7i,72,73,0) - — J7,-i(7i, 72, 73, 0)] 
71 71 

-^^(7? + 72 + 73')[ C^jKti, 72,73, - -C/fi(7i,72,73,??)] + 

71 

(3.14) +[i^|(7i, 72,73, r/) - -F,!(7i,72,73,ry)] 

71 

'7[ C^,|(7i,72,73,?7) - —[^n (71, 72, 73,??)] - [ C7j-3(7i, 72, 73, 0) - — C/,-i(7i, 72, 73, 0)] 
71 71 

-^7? + 72 + 73')[ £^,1(71,72,73,??) - -C/fi(7i,72,73,??)] + 

71 

(3.15) + [ ^;K7i,72,73,??) - -^;!(7i,72,73,?/)] 

71 



(3.16) 7i[7j^(7i,72,73,f?) + 72 [7^(71,72,73, 77) + 73 [7,® (71 , 72 , 73 , '?) = 

(3.17) [/ji(7i,72,73,0) = [71(71,72,73) 

(3.18) [^,2(71,72,73,0) = [72°(7i,72,73) 

(3.19) [/j-3(7i,72,73,0) = [73(71,72,73) 

Let us rewrite system of equations p.l4p — ()3.16p in the following form: 
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(3.20) 



[ —^i! (71, 72,73,'?) - -^"1(71, 72, 73,??)] +[ — C/,i(7i,72,73,0) - t/,-2(7i, 72, 73, 0)] 
71 71 



77 + i^ili + li + li)] — C/,® (71, 72,73, T?) - [v + J^(7? + 7I + 7I)] C/|(7i, 72,73,??) = 



71 



(3.21) [ ^i^j^ (71, 72,73, '/) - ^j|(7i,72,73,??)] + [ ^C/,i(7i,72,73,0) - C/,-3(7i, 72, 73, 0)] 



(3.22) 



7it^j^(7i,72,73,??) + 72 ?^j|(7i,72,73,»7) + 73 C^j|(7i, 72, 73, = 



Determinant of this system is 

' V + + 7I + 7I)] + + ^2 + 7I)] 



A 



[ ^ + i'(7?+7l+7l)] ^ 
71 





72 





^ + K7?+7|+7|)] 
73 



(3.23) 



T] + v(nl + 7I + 7I)] (7? + 7I + 7I) 

71 



^ 



And consequently the system of equations p.l4p — (|3.16p and/or p.20p — (|3.22p has a unique 
sohition. Taking into account formulas p.l7p — (|3.19p wc can write this solution in the following form: 



t^ji(7i,72,73,??) 



(3.24) 



[(7I + iVlFfx (71 , 72 , 73 , ^7) - 7i72 (71 , 72 , 73 , ^7) - 7173^^1 (71 , 72 , 73 , ??)] 
(7? + ll + 7i) + + 7I + 

C/i(7i,72,73) 



N + i^(7f + 7I+7I)] 



t^j|(7i, 72,73,??) 



(3.25) 



[(73 + 7i (71 , 72 , 73 , ??) - lilzFf^ (71 , 72 , 73 , ??) - 727i J^j^ (71 , 72 , 73 , ??)] 
(7? + ll + 7I) [?? + + 7I + 7I)] 

1^2(71,72,73) 



+ 



t{,l(7i,72,73,??) 



(3.26) 



[77 + i^(7? + 7|+7|)] 

[(7i +72)-Pj3(7i,72,73,??) -737i-Pji(7i,72,73,??) -7372-^^1(71,72,73,??)] 
(7? + ll + 7I) [?? + + 7I + 7I)] 

1^3(71,72,73) 



+ 



[?? + i^(7? + 7|+7|)] 
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Then we use the convolution theorem with the convolution formula (jA.4[) and integral (jA.5[) for p.24[) 
p.26p and obtain: 



(3.27) 



C/,-i(7i,72,73,i) = 

,(^j'+-y;+-yj)(f-r) [(7I +7l)-F3i(7i, 72,73, t) - 7i72fj2(7i,72,73,T) - 7173-^,3(71: 72, 73, t)] 

(7? + 7I + 7I) 

+ e -''(^i'+^^'+^^'^*[/?(7i,72,73) 



dT 



C/?-2(7l,72,73,t) 



[(7I +7i)-P)2(7i,72,73,t) - 7273-Fj3 (71, 72,73, t) - 7271-^)1(71, 72, 73, t)] 



(3.28) 



ill + ll + 7|) 
+ e -''(^^'+■^^'+^3')* [7° (71, 72, 73) 



dT 



C7j3(7i,72,73,i) = 

g _^(^2^,^2^^2^(^_^) [(7i +7l)-Fj3(7i, 72,73, t) -737i^ji(7i,72,73,T) - 7372^32(71, 72, 73, t)] 
ill + 7I + 7|) 



(3.29) 



+ e -(7, +7. +73)* [/0(7i,72,73) 



-Pi(7i, 72, 73, is obtained from equations (|3.ip . p.2p , (|3.3p with use of equation p.4p 



(3.30) P,(7i,72,73,0 = i 



[7i-Fji(7i,72,73,0 +72-Fj2(7i,72,73,^) + 73-Fj3(7i, 72, 73, ^)] 
(7?+7l+7|) 



Using of the Fourier inversion formula (jA.ip leads to: 



Uji(a;i,X2,X3,i) 



(27r)3/2 7_ 



00 /'OO t^OO 



00 J —00 J —00 



g -,.(y2^y2^+y2)(t-r) [(7l + 7g)Jjl(7l,72,73,T)] 

(7? + 7l+7l) 



e -'^(7?+7?+7f)(t-r) [7172-^32(71, 72, 73, t) + 7i73-Fj-3(7i, 72, 73, t)] _^ 

(7?+7l+73) 



+ e -''W+^=+^3)*f/°(7i,72,73) 



1 

8^ 



00 poo poc / 2 I _ 2 



00 ^ — 00 J —00 



(.li + 73) 

(7?+7l+7|) 



g -''(7r+72+73)(*— ^) 



g -»(^l7l+^272+^373) (^^^^^2^73 = 

g •i(:ri7i+:r272+i373) 



00 pOO pOO 



00 ^ —00 ^ —00 
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■fjl{xi,X2,X3,T)dxidx2dx3dT 



-i(£Ci 71 +2:272 +2:373) 



1 

8^ 



1 

8^ 



00 />00 />00 



00 —00 J —00 



7172 



(7? + 7I + 7I) 



e -i^(7r+75+73)(t-^) 



00 /'CXD />00 



00 ^ — 00 ^ —00 



^71^72^73 - 
g 1(2171+2272+5373) 



■/j-2(£i,i2,i3,''') dxidx2dx3dT 



g -,(2:171+2:272+2:373) d-fid-i2dl3 



00 /'OO /'OO 



— 00 ^ — 00 



7173 



g -''(71' +72 +73 )(*-^) 



-00 (7? + 7I + 7|) 

■fjz{ii-,i2,X3,T) dxidx2dx3dT 



00 /"CXD />00 



— 00 ^—00 J —00 



1(2:171+2272+5373) 



g -,(2:171+2272+2:373) d-iid-i2di3 + 



1 

8^ 



OO /'OO fOC 



00 J ~oo J —00 



g -"^(71+72+75)* 



00 /"OO /"CJO 



00 J —00 J ~oo 



g 4(2171+5272+2:373) . 



Ul{xi,X2, X3) dxidX2dX3 



g -,(2171+2272+2:373) d-fid-f2d-f3 



(3.31) 



Uj2{xi,X2,X3,t) 



(27r)3/2 7_ 



00 poo poo 



00 J —00 J — 00 



^(7i+72+73)(*-^) 

[(7|+7?)^J2(71 ,72,73,TjJ 

(7? + 7I + 7I) 



_ /"*g -,.(7?+7,"+7.?)(t-r) [7273^33(71: 72, 73, t) + 7271 -P'jl (7I: 72, 73, t)] 



+ e -''(^i'+^^'+^3)*[/°(7i,72,73) 



(7?+7l+7|) 

g -4(2171+2272+2:373) dj^dj2dj3 



1 

8^ 



00 poo poo 



00 'J —00 'J —00 



7271 



(7?+7l+7|) 



g -''(7i:+72+73)(*— ^) 



00 /'OO roo 



— 00 —00 —00 



4(2171+2272+5373) 



■fjl{xi,X2,X3,T) dxidx2dx3dT 



g -,(xi7i+2272+2373) djidj2dj3 + 



87r3 



1 

8^ 



1 (^32+^2) 



— 00 ^ — 00 



00 poo pOO 



g -i^(7?+72+7|)(*-T) 



— 00 ^ — 00 



-00 (7?+7|+7|) 

•/j2(ii, 52,^3,''') dxidx2dx3dT 

^273 r /"^g -,(^^3+^^2+^|)(t-,) 

-00 (7?+7|+7|) Uo 

•/j3(il,X2,i3,T) dxidx2dx3dT 



00 /'OO /'OO 



— 00 —00 —00 



4(2171+2:272+5373) 



g -,(2:171+2272+2373) djidj2dj3 



00 /'OO /'OO 



— 00 —00 —00 



4(2171+2272+5373) 



-4(2:171+3:272+2373) 



d'yid'y2dj3 
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1 

8^ 



CXj J ~00 J — OO 



OO poo pOQ 



OO J ~oo J ~oo 



i(5l7l+£272+5373) 



W2(ii, a;2, is) dxidx2dxz 



(3.32) 



52l(/,l) + 522(/,2) + 523(/,3) + B{ul) 



(27r)3/2 



OO nOO poo 



OO J — OO J ~oo 



e -'^(7?+7.^+7g)(t-r) [(7? + 7|)-Fj3(7i,72,73,t)] 

(7? + 7l+7l) 



e -i^(7r+72+73)(t-^) 



I [737i-Fii(7i,72,73,T) +7372i^j2(7i,72,73,T)] 



dr + 



+ e ^''(^'+^'+^')*f/3°(7l,72,73) 



(7?+7l+7i) 

g -^(^171+^272+^373) ^^^^^2^73 



1 

8^3 /_ 



OO /'OO /'OO 



OO ^ — OO 



7371 



g -''(7r+72+73)(*-^) 



-OO {ll + ll + 7l) 

■/ji(£i,i2,53,T) dx\dx2dx-idT 



OO /"OO />oo 



— OO OO J — OO 



i(5 171+^2 72 +5373) 



-1(3:171+3:272+3:373) 



(3.33) 



1 

8^ 



1 

8^ 



00 /'OO /'OO 



7372 



00 J —00 ^ —00 



00 /'OO /'OO 



(7f +71+73) 



g -i^(7r+72+73)(*-'^) 



00 poo /'OO 



00 ^ —00 —00 



^71^^72^73 - 

g i(5l7l+i272+5373) 



■/j2(ii,i2,i3,T) dxidx2dx^dT 



00 J ~oo J —00 



(7?+7l) 
(7? + 7I + 7,1) 



g -''(7r+72+7j)(t-T) 



g -^(3:171+3:272+3:373) d'yid'y2d-/3 4 

g i(ii7i+S272+i373) 



00 /'OO /"OO 



00 J —00 J —00 



■fj3ixi,X2,X3,T) dxidx2dx3dT 



g -^(3:l7l+3:272+3:373) ^^^^^2^73 + 



1 

8^ 



00 pOO pOO 



00 ^ — 00 ^ —00 



g -"^(7r+72+73)* 



00 poo poo 



00 J —00 J ~oo 



g i(Sl7l +5272+^373) . 



u^{xi, X2, X3) dxidx2dx3 



g -,(3:171+3:272+3:373) ^^^^^2^73 



53l(/,l) + 532(/,2) + 533(/,3) + B{ul) 



Pj {xi,X2,X3,t) 



(27r)3/2 



00 poo poo 



00 ^ —00 ^ —00 



[71 Fj 1 (71 : 72 , 73 , ^) + 72 -Fj 2 (71 , 72 , 73 , 
(7? + I2 + 7|) 
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8^ 



I 

8^ 



I 

8^ 



oo /"oo poo 



oo ■/ — oo ■/ — oo 



(7?+7l+7|) . 
7i 

(7? + 7| + 7|) 



00 /"OO /"OO 



■fji{xi,X2,X3,t)dxidx2dx3 

72 

(7? + 7I + 7I) 



00 /'OO /'OO 



00 J —00 ^ — 00 



00 —00 •/ —00 



00 /'OO nOO 



• fj2 (xi , X2 , , i ) dxidx2 dx^ 

73 

(7? + 7| + 7|) 



00 /'OO /'OO 



00 ^ — 00 ^ — 00 



(3.34) 



■/33(ii,i2,^3,0 (ixi(ii2(i^3 



5l(/,l) + S2{f,2) + 53(/,3) 



-^(:ri7l+:^272+:r373) djidj2dj3 = 
g i(Sl7l+5272+S373) 

g -»(a,i7l+^272+^373)^^^^^2^^g ^_ 

g i(il7l+X272+S373) 

g -^(:.i7i+a,272+:r373) ^^^^^2^73 + 
g i(ii7i +£272+^373) . 

g -i(xi7i+.,272+:r373) ^^^(^^2(^73 = 



00 •/ —00 •/ —00 



00 /'OO /'OO 



00 — 00 ^ —00 



So, the integrals (j3.3ip — (|3.34p exist [look below in a), b) - the Fourier transform of the elass S, 
infinitely differentiable functions]. For these caleulations inverse Fourier transforms are defined as Cauchy 
principal values and for 7i=0, 72=0,73=0. 

Here 5n(), 5i2(), ^laO, 52i(), ^22 0, ^23(), ^aiO, ^320, ^aaO, B(), 5i(), 52(), SaO are 

the integral - operators. 

5l2() = 521 

5i3() = 5ai() 

523 - 5a2() 

We have for the vector Uj from the equations (|3.31l) — (|3.33l) : 



(3.35) Uj = S ■ fj + B - if , 

where S and B are the matrix integral operators: 



5ii 


5l2 


5i3 


^21 


S22 


523 


531 


S32 


5a3 



/BO \ 

B 
V OBJ 
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and for the funcion pj from the equation p.34p ; 

(3.36) p, = ~S ■ f,, 

where S is the matrix integral operator: 



Si 






















^3 



We put fj from equation (j2.20p into equation (|3.35p and have: 

= 5 ■ ( / - ( • \7)uj-i) + S-u" = 

= S- f - B- (uj_i • V ) Uj-i + B-iP = 

(3.37) ^ ui — S ■ {uj-i ■ V) Uj-i 

Here ui is the solution of the system of equations (|2.13p — (|2.20p with condition: 

3 



E""?^ = k=l,2,3 



n=l 

For j = 1 formula p.35p can be written as follows: 

(3.38) ui ^ S ■ h + B-iP , h{x,t) - f{x,t) 

If t — 5- then ui ^ vP (look at matrix integral operators S, BQ - integrals p.3ip — (|3.33p ). 
For j = 2 we define from equation p.20p : 

(3.39) f2{x,t) = fiix,t) - {ui ■ V)ui 
We denote: 

(3.40) f; = {ui ■ V) ui 
and then we have: 

(3.41) f2ix,t) = fi{x,t) ~f; 
Then we get U2 from (|3:35l) . ([XSg]) : 

(3.42) U2 = S ■ f2 + B -if = S ■ ifi ~ f;) + B -if = ui - U2 



NEW APPROACH. SOLUTION AND ITS JUSTIFICATION. 



Here we have: 



(3.43) ul = S ■ /* 



If t — ^ then ^ (look at matrix integral operator 5* - integrals p.3ip — p.33|) ). 
Continue for j = 3. We define from equation p.20|) : 



(3.44) /3(x,t) = /i(x,t) - {U2 ■V)u2 
Here we have: 

(3.45) • V) U2 = ((ui ~ ni) {u, - ui) = f; + f; 

We denote in (|3.45p : 



(3.46) /a = - (ui ■ V) U2 - (u^ • V) mi + (m^ • V) u\ 
and then we have: 

(3.47) Mx,t) = /i(x,t) -/t - it 
Then we get U3 from ([SlSS]) . ([S^S]) . ([OH]) : 

(3.48) = § • /I + S • ^0 = § • (/i - - it) + ^ • = "1 - "2 - " 
Here we denote: 

(3.49) - § • it 

If t — > then -1*3 — !■ (look at matrix integral operator S - integrals p.3ip — (13.33^ ). 
For j = 4. We define from equation (|2.20p : 

(3.50) U{x,t) = hix,t) - (7X3 • V)m3 
Here we have: 

(3.51) (7X3 • V) ^^3 = {{U2 - W3) • V ) - ui) = f; + f; + n 

We denote in (P3T|) : 

(3.52) f* — ^ {u2 ■ V) U3 - [u*^ • V) M2 + (M3 ■ V) U3 
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and then we have: 

(3.53) Ux,t) = /i(.T,i) -f; -f* -fl 
Then we get U4 from ((3:35)) . ([XSH)) . ([OH)) . dSH)): 

(3.54) ^4 = 5 • (/i - it - it - A*) + B-^^ = ^ii - ul - ul - ul 
Here wc denote: 

(3.55) = 5 . /* 

If t ^> then ■U4 ^ (look at matrix integral operator S - integrals p.3ip — ()3.33p '). 
For arbitrary number j {j >2). We define from equation ()2.20p : 

(3.56) fjix,t) = hix,t) - iu,.i ■ V)w,_i 
Here we have: 



(3.57) {u,^, ■ V)u,_, =T.f* 

and it follows: 



1=2 



(3.58) f, ^ h ~Y.fl 



1=2 



Here /| is taken from formula p.40p and 

(3.59) ./t = - ■ V) uU - ("f-i • V) + {ut-i ■ V) uU (l > 2) 
Then we get uj from (|335)) . (jXSSl 

j 

(3.60) = § • i- + B-^ZO = 5 • (/i - ^/T) + 5-72" = ^Zi - 

1=2 1=2 

Here we denote: 

(3.61) u*i =~S ■ f: {2 < I < j) 

If t — !> then Ul ^ (look at matrix integral operator S - integrals ()3.3ip — ()3.33p '). 
We consider the equations p.38|) - p.6ip and see that the series p.60p converges for j -> 00 
with the conditions for the first step (j = 1) of the iterative process: 
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Y^Uo^^^O k = l,2,3 

n—1 

Hence, we receive from equation p.37p when j ^ oo: 

(3.62) Uoo ^ Ui - S ■ {Uoo ■ V) Uoo 

Equation (|3.62p describes the converging iterative process. 
Below we show a proof that the iterative process is converging. 

a) Let us consider S to be the class of all infinitely differentiable functions (p{x) {~oo < x < oo), 
satisfying inequalities of the form 

(3.63) Ix'^ip'^'^^x) \< Cfcgforanyk,q = 0,l,2,... 
where Ckq is a constant and depends on f{x). 

Then, FS = S, i.e., the Fourier transform operator F maps the class S onto the whole class S [TB] 
Now let us rewrite conditions (|2.17p , (|2.18p in the following form: 

(3.64) I (1+ I i \)^d^v°{i) I < Cc^K on for any a , K 



(3.65) I (1+ I i I +T)^9?9f /(i, r) | < C^pK on x [0, oo) for any a , ^ , K 
or for arbitrary k (1 < fc < N) 

(3.66) I (1+ I i D^d^uUx) I < CaK on R^ for any a , K 

(3.67) I (1+ I ^ I +T)^a?a^/fc(i, r) I < Ca/3K on R^ x [0, oo) for any a , ^ , K 

By comparing p.66p , p.67p with p.63p we can see that infinitely differentiable functions u°(i), /fc(S, r) 
are satisfying inequalities of the type (|3.63p and hence, w°(i) £ S{R^), /fc(i,r) G S{R^). 

Let us consider the first step of iterative process. 

We can see that inner integrals (Fourier transforms) in the integral operators SnQ, S12O, SisQ, S21O, 
522(),523(),53i(),532(),^33(),S() from formulas (0311) , (IS , (E331) transform ^0(5), /^(x, r) into 
Mg(7) e SiR^), /fe(7,T) e S{R^), according to [16]. 

Muhiplication oiuHj) by 6 -"^(ti'+t'+^D* and /fc(7,r) by 6 -'^(t?+72+7|)(*-^) and by fractions 
I (T2+7I) 1^ 1 I (71-72) 1^ 1 I (71-73) 1^ 1 

I (7^+72^+73^) l< ^' I (7^+72^+71) l< ^' I (7?+7|+7|) ^' 
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I (T3+7?) 1^ 1 I {72-73) 1^ 1 I (71+72) 1^ 1 
I (7^+72^+73^) l< ^' I (7^+7|+7|) l< ^' I hi+7i+7i) l< ^ 

keeps result functions in class S{R^). 

Inverse Fourier transforms (outer integrals in the integral-operators 5'ii(), 5120, ■S'laO, S'2i(), 
S'22(),523(), 5*310, 532(),533(),B()) result in uik{x,T) e S{R^) according to [l6]. Integrating uik{x,T) 
with respect to r over the interval [0, t] keeps functions uik{x,t) in class S{R^). 

Let us consider the second step of iterative process. 
We obtain the first correction to applied force 

(3.68) it = {ui ■ V) ui 
from formula (|3.40p . For arbitrary k (1 < fc < N) we have 

(3.69) ^ yu, 



^1 



From uik{x,T) G 5(i?^) it follows that ^ G 5(i?^), and hence 7*^ G 5(E^) [H]. 
We can obtain the first correction to velocity 

(3.70) u* = ~S ■ 

from formula (j3.43p . After analogous reasoning for components Ujj, like for uik on the first step of 
iterative process, we have u^/. G S{R^) according to [16]. 

Hence, we have received that on any arbitrary step / (/ > 1) of the iterative process a correction to the 
force /;* as well as correction to the velocity W;* are infinitely differentiable functions and fj^). G S{R^), 
G 5(i?^) (1 < fc < TV). 

b) Following [TB] we introduce classes of functions Wm and in this paragraph. Let M(x) and ^}{t) 
be dual functions, in Young's sense, and let Wm be the class of all infinitely differentiable functions 
f{x) (—00 < a; < 00), satisfying the inequalities 

(3.71) |^(«)(x) I < C,e (g = 0,l,2,...). 

where is a constant and depends on (p{x). 
If iIj{s) = F[ip{x)] is Fourier transform, then 

(3.72) I + *r) I < C'^e {q = 0, 1, 2, ...). 

Let W^^ be the class of all entire functions V'(s) satisfying inequalities of the form ()3.72p . 

Then, FWm = W^, in other words, the Fourier transform operator F maps the class Wm onto the class 

W^ and FW^^ = Wm, i-C-, the Fourier transform operator F maps the class W^ onto the class Wm jl6) . 
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Now let us consider u^{x) S Wm{R^), fk{S:,T) £ Wm{R'^) and go to the first step of iterative process. 
Inner integrals (Fourier transforms) in the integral-operators S'ii(), 5120, SisQ, S'2i(), 'S'22(), 'S'23(), S31O, 
S'32(),S'33(),S() from formulas (la^Tll , (l332l) , (jSlall transform ^^(i), /fe(i, r) into ^^(7) g W^\R^), 
fkh,T) e W^\R'^), according to [H]. 

Multiplication of ^0.(7) by e -'^(T?+72+7|)t and 7^(7, r) by e -''(7?+72+73)(t-r) ^^d by fractions 
I (72+73) 1/ 1 I (71 ■72) 1^ -1 I (71 ■73) 1^ 1 

I (^?+^|+^ l< 1' I (7?+72^+7|) l< 1' I (7?+72^+7|) l< 1' 

(73+7?) 1^ 1 I (72-73) 1^ 1 I (7?+72) 1^ 1 
(71^+72^ +73') l< ^' I (^?+^^?+^ l< ^' I (71^+72^ +73') l< 

keeps result functions in class W^{R'^). 

Inverse Fourier transforms (outer integrals in the integral-operators Sii{), S12O, S13O, '5'2i(), 

8220^ S23O, S31O, S32O, SssO, B{)) result in Mifc(a;,T) e WuiR'^) according to [TS|. Integrating mu- (a;, r) 

with respect to r over the interval [0, t] keeps functions uik{x,t) in class Wm{R^)- 

Let us consider the second step of iterative process. 
We obtain the first correction to applied force 

(3.73) f; = (ui ■ V) ui 
from formula (|3.40p . For arbitrary k (1 < fc < N) we have 

(3.74) = E--?^ 

From uik{x,T) e Wm{R'^) it follows that ^ G WmIR^), and hence f;,^ G W^M(i?^) [S]. 
We can obtain the first correction to velocity 

(3.75) ul = S ■ 

from formula (j3.43p . After analogous reasoning for components Ujj, like for uik on the first step of 
iterative process, we have Ujfc ^ Wm{R^) according to [16]. 

Hence, we have received that on any arbitrary step I {I > 1) of the iterative process a correction to the 
force /;* as well as correction to the velocity M;* are infinitely differentiable functions and G Wm{R^), 



u 



Ik 



G Wm{R^) (1 < fc < N). 



c) Let us estimate superiorly solution of the Cauchy problem for the 3D Navier - Stokes equations 
by iterative method. The purposes of this estimation are: 

1) to show convergence of the iterative method; 

2) to obtain analytical form of the first and second steps of the iterative process; 

3) to receive estimated formula for the border of convergence region of the iterative process in the 
space of system parameters. 

We substitute fractions 
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(7I + 7I) I -, I (71 • 72) I , I (71 ■ 73) 

(7? + 7| + 7|) ' ' ' (7? + 7| + 7|) ' ' ' (7? + 7| + 7|) 



^ ^ ^ ' (7? + 7l + 7l) ' ' ' (7? + 7l + 7l) ' ' ' (7? + 7l + 7l) ' 

by 1 in t he in te gral-o perators S'nO, S'i2(), 5*130, 5210, 'S'22(), 'S'23(), S'aiO, 5*320, S'saO from formulas 
(|3.3ip . (|3.32p . (|3.33p for all steps of iterative process. 
Then we take 

/ii(ii,i2,i3,T) = F-/(r)-e -m'(^?+S2+^1)^ F, fi - constants, F > 0, fi > 

(3.77) /12 = /i3 = 0, iTo = 

After Fourier transforms (inner integrals in the integral-operators ^nO, 5*210, '^'siO from formulas 
((33T1) . ((3321) . dSSl) we have: 

/ \ 3/2 (-j,2_^_^2^^2 j 

(3.78) /ii(7i,72,73,r) = ^ ' /(^) ' f " ' ^'"^ ' 

Now wc multiply /ii(7i, 72, 73, t) by 6 ~''(ti+^2+T3)(*-'^)^ change order of integration by t and 71, 72, 73 
and after Inverse Fourier transforms (outer integrals in the integral-operators 5ii(), 52i(), 53i()) we have: 

F ■ fir) -^"(x^+x^+xg) 
^n(.i,X2,.3,r) = [,^2,(,_;);i]3/2 -e - > 

(3.79) Mil = -"12, "12 = "13 

Then we get: 



/(r) 



(3.80) ?2ii(xi,X2,X3,i) = [V^(/^V) + l]3/2 -^ ^I^lJIT^dr 

We substitute y for r: y = [4^2^(t-T)+i] ' [4p^4t'-^)-H]^ •^'^^^ ^/^^^ and receive after 

integration: 



?2ii(xi,X2,X3,t) ^ ^-^-^-^—^—^ 



2, 2 2 2^ '-1+-2+-3' 

_ g -A* (2:1-1-2:2+3:3) _|_ g (4^:57:7+1) 



Aii-^v(x\ +x1+ x|) 



(3.81) wii = -"12, "12 = "13 
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Here j{a, x) is the incomplete gamma function |31j and $(a, c; x) is a confluent hypergeometric function 



m 



duii{xi,X2,X3,t) 
dXn 



F ■ Xr. 



1 



2/x4jy(a;2 + x| + xl) \ {xl + xl+xl 

1 



2 / 2 I 2 I 2\ 



2 / 2 I 2 I 2\ 

g -A* (a:i+a:2+a;3) 



(4/i2zyt + 1) 



g (4m^i'* + 1) 



(3.82) 



F • Xfi 



^2,i■~^l\xi + xi+xi)] 



{A^i^iyt + 1) 



' ni-*[2,3;^4i^±tf|)]^ (l<n<Ar) 



(4^2 _^ 1) 



Here ui{xi, X2, x^^t) and '^^'^i^y-^^-^^^-'^) g^j-g estimations from above of velocity and partial derivatives of 
velocity on the first step of the iterative process. 

Now we do a next level of superior estimation of velocity and partial derivatives of velocity to continue 
the iterative process. For this goal we receive from formulas p.8ip . p.82p with condition {x\ +x\ + x^) > 1 



(3.83) 



Uii{xi,X2,X3,t) 



e 



(3.84) 



duii{xi,X2,x^,t) F 



' dxn 2^i'^v {A^j?vt + 1) 

Then in this case we have from formula p.74p following estimation from above: 



(3.85) 



I /2I I = I /22 I - I /23 I = I "11 



duii 



F2 



1 



-2,.^(x2 + x^ + x^) 



dXn 4/i4l/2 {A^l'^vt+l) 



After Fourier transforms (inner integrals in the integral-operators 5ii(), 5*210, S^ii) from formulas 
dSSU, ((332|) . dSjSl)) we have: 



/2fe(7l,72,73,'r) 



F' 



1 



7r(4^2j^^ _l_ 



(3.86) 



(1 < fc < A^) 



3/2 



• e 



(4f.^,^x + l)(7^+T.|+7|) 



Now we multiply /2*fc(7i, 72, 73, t) by 6 -''(Tr+T2+73)(*-^), apply substitution of fractions ((3J6)) by 1, 
change order of integration by r and 71 , 72 , 73 and after Inverse Fourier transforms (outer integrals in the 
integral-operators 5'ii(), S'2i(), 5*310) we have: 



'«2fc(a;i,a;2,X3,T) 



F2(4/i2l.T+ 1)1/2 



4/i4l/2[8/i2;y(t - r) + (4/i2jyr + l)]3/2 



(3.87) 



(1 < yfc < iV) 
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Let us further increase level of estimation. To do so we substitute G [8M2''(t-T)+(4M2^T+i)] from l|3.87p by 



Then we get: 

(3.88) ii;.(...„X3,.) = — - (8M-l'-vJ+l)3/^ -^- 

We substitute y for r: y = (s^Vt-l^iVr+i) ' = {sTI^l^t^f^^W+TF^'^ ^^"^ receive after integration: 



(3.89) u;^ixi,X2,X3,t) = ^^-e (^^^-+1) 



^2 -M^(xf+x^+x^) 

8/A.3 



1 - 7r/4 - — , + arctg- 



(8/ii/t+ 1)1/2 ^(8^1/^+1)1/2 



Now wc compare uii{xi, X2, X3,t) from formula p.83p with W2j,(xi,X2,X3,t) from formula p.89p and 
see that the iterative process is converging with estimated condition: 

(3.90) ^ < 1 

where F and /i were introduced in formula (|3.77p . Condition (|3.90p is the estimated formula for the 
border of convergence region of the iterative process in the space of system parameters. 

For arbitrary step j (j > 2) of iterative process we may take Uj from formula p.60p and apply estimation 
algorithm analogous to formulas (|3.76p - p.89p . 

Since these results are shown for the superior estimation, then for precise calculations the convergence 
of the iterative method will be even better. 

Then we have from formula (|3.34p : 

(3.91) Poo - SiiUl) + ^2(/oo2) + ^3(/oo3) 

Here foe = (/ooi, /oo2, /cxjs) is received from formula p.58p . 

On the other hand we can transform the original system of differential equations (j2.7p — p.9p to the 
equivalent system of integral equations by the scheme of iterative process p.35p , p.37p for vector u: 

(3.92) u = ui ^ §-{u ■ V) u, 

where ui is from formula p.38p . We compare the equations p.62p and p.92p and see that the iterative 
process ()3.62p converges to the solution of the system p.92p and hence to the solution of the differential 
equations dSJ]) - ([23]) . 

In other words there exist smooth functions Poo(x, t), Uooi(x, t) (i = 1, 2, 3) on x [0,00) 
that satisfy dH]), (H^), and 



(3.93) Poo, Uooi e C°°(R3 X [0, 00)), 
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(3.94) 

for all t > 0. 



/ |Uo„(x,t)pdx<C 

JR3 



4. Justification of the analytical iterative method solution for Cauchy problem for the 3D 
Navier-Stokes equations. 

4.1. Spaces S, TSt [20], [27]. 

Let us consider space S of all infinitely differentiable functions </?(x) defined in N-diniensional space 
(N = 3), such that when | a; | — >■ oo these functions tend to 0, as well as their derivatives of any order, 
more rapidly than any power of j^j. 

To define topology in the space S let us introduce countable system of norms 



(4.1) II^IIj, = sup, \x''D'>if{x)\ (p- 0,1,2,...) 

\k\, \q\ <p 

where 

\x'^DM-)\ = \4'----%"'-^^^^r^S^\ 

fc = (fci, . . . ,fcAr), (7= (qi,...,(7Ar), x'' = x'^^ . . . x']:^ 
= olfJZ- ' (fci,...,<7w = 0,l,2,...) 

Space S is a perfect space (complete countably normed space, in which the bounded sets are compact). 
Space of vector- functions (f is a. direct sum of N perfect spaces S (N = 3) [H], [15] : 

To define topology in the space let us introduce countable system of norms 



N N 

(4.2) = Y,y:\\p '""P- l^'^V.(^)l (p = 0,l,2,...), (iV = 3) 

1=1 ,= i\k\.k\<P 

4.2. Equivalence of Cauchy problem in differential form (j2.ip - (|2.3p and in the form of an 
integral equation. 

Let us denote solution of the problem p.ip - p.3p as {u{xi,X2,X3,t), F{xi,X2,X3, t)}, in other words 
let us consider infinitely differentiable by t G [0,oo) vector-function u{xi, X2, X3,t) G and infinitely 
differentiable function P(a;i, a;2, X3, t) G S, that turn equations p.ip , ()2.2p into identities. Vector-function 
u{xi,X2,X3,t) also satisfies the initial condition p.3p {if'{xi,X2,X3) G T^): 



(4.3) U{xi, X2,X3, t)\t=n = {l°{xi,X2,X3) 

Let us put {u{xi,X2,X3^t), P(xi, 2:2, 2^3, t)} in to equations (j2.ip , p.2|) and apply Fourier and Laplace 
transforms to the result identities considering initial condition ()2.3|) . After all required operations (as in 
parts 2 and 3) we receive that vector- function u(xi,X2,X3,t) satisfies integral equation: 
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(4.4) u ^ -S ■{u-V)u + ui = S -u 
where uieT^ is from (jX^ : 

(4.5) Ml = S-f + B-v° , 

Here / G T§, vP G T\^. Results of operators S ■ B ■ if, S ■ {u- V)it are also belong T~§ since Fourier 
transform maps perfect space f:^ onto f:^. Function P G S is defined by formula 

(4.6) P = S-f- S-{u-V)u, 

where vector- function u is received from (|4.4p . 
Here S, S , B, S are matrix integral operators. 

Going from the other side, let us assume that u{xi,X2,X3,t) G is continuous in t G [0,oo) solution of 
integral equation (|4.4p . Integral-operators Sij • [u • V)u are continuous in t G [0,(X)). From here we receive 
that according to (jHH), (|4?5|) 

u(a;i,X2,a;3,0) = vP{xi,X2,xz) 

and also that u{xi,X2,x^,t) is differentiable by t G [0,cxd). As described before, the Fourier transform 
maps perfect space on itself. Hence, 2:2, X3, t) and P{xi, X2, X3, t) from formula (|4.6p is the solution 
of the Cauchy problem (|2.1I) - (|2.3|) . From here we see that solving the Cauchy problem (j2.ip - (|2.3p is 
equivalent to finding continuous in t G [O.cx)) solution of integral equation ()4.4p with condition ()4.5p . 

4.3. The fixed point principle. [21] , [22], [23], [24], [25], [26] 

Let us use the fixed point principle to prove existence and uniqueness of the solution of integral equation 

631). 

For this purpose we will operate with the following properties of matrix integral operator S : 

1. Matrix integral operator S continuously depends on its parameter t G [0,oo) (based on formulas 
(l?3TD - dSSSD). 

2. Matrix integral operator S maps vector-functions u from perfect space onto perfect space T^. 
This property directly follows from the properties of Fourier transform [20| , and the form of integrands of 
integral operators Stj, B (based on formulas (|3.3ip - (|3.33p and paragraphs a) and b) in part 3). 

3. Matrix integral operator S is "quadratic" and therefore we will consider vector- functions ui, u such 
that I Ml |< 1, I M |< 1 for any value t G [0,oo). 

For example, the convergence of the iterative process in part 3 is achieved because of condition p.90p : 



4. US' ■ u ~ S ■ u \\p < \\u — u lip for any u, u £ TS {u ^ u , use property 3) and any t G [0,c») 
(based on formulas p.3ip - (|3.33p ). 

Let us prove that matrix integral operator S is a contraction operator. 
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Theorem 1. Contraction operator. 

Matrix integral operator S maps perfect space onto perfect space and for any u, u (z 
(u ^ u ) satisfying property 3, the condition 4 is valid. 

Then matrix integral operator S is a contraction operator, i.e. the following condition is true: 

(4.7) US' ■ u ~ S ■ u \\p < a ■ \\u — u \\p 
where a < 1 and is independent from it, u G for any t G [0,oo). 
Proof by contradiction. 

Let us assume that the opposite is true. Then there exist such Un, G (n=l,2,. . .), that 

(4.8) \\S -Un - S -uJlp = a„ • - uj\p (n = l,2,...; a„ ^ 1) 

Also, because is a perfect space, we can consider that u„ — > u G and ^ u E T^. Then 
the limiting result in (|4.8I) would lead to equality 

US' -u - S ■ u \\p = \\u - u lip , 

which is contradicting condition 4. Hence, S is a contraction operator. 

Theorem 2. Existence and uniqueness of solution. j21j 

= V . . . ^ . r-7— n, 

Let us consider a contraction operator S . Then there exists a unique solution u* of equation (|4.4p in 
space for any t G [0,cxj). Also in this case it is possible to obtain u* as a limit of sequence {u„} , 
where 

Un+l = S {un) (n = 0, 1, . . .), 
and uq = 0, and ui is taken from (|4.5p . 

The rate of conversion of the sequence {un} to the solution can be defined from the following inequality: 

a" 

(4.9) ||u„ - u*\\p < -\\ui - uoWp (n = 0,l,...) 

(1 - a) 



Proof: 
Since 



then according to (|4.7p : 



Un+l = S (u„), Un = S (u„_i) 



llWri+l - Un\\p < a-\\Un - Un-lUp- 

Subsequently using analogous inequalities while decreasing n we will receive: 

||m„+1 - Un\\p < a" ■ ||wi - WoIIp- 
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From this result it follows that 



\\Un+l - UnWp < \\Un+l - U„+(_i||pH h||?i„ + l - U^Wp < 

it" 

(4.10) < (a"+'-i + • • • + a")||ui - uo||p < jz r\\ui - uoWp- 

(1 - a) 

Because of a" -> for n — oo, the obtained estimation ()4.10p shows that sequence {w„} is a Cauchy 
sequence, and since the space is a perfect space, this sequence converges to an element u* 
such that 5* {u*) makes sense. We use inequality (|4.7p again and have: 

\\un+i - S^{u*)\\p = WS^iun) - S^iu*)\\p < a-\\un - U*\\p (n = 0, 1, 2, . . .) 

The right part of the above inequality tends to for n — > cxd. and it means that Un+i S (u*) 

and u* ^ S (u*)- In other words u* is the solution of equation (|4.4p . 

Uniqueness of the solution also follows from (|4.7p . In fact, if there would exist another solution 

then 

|l5 - u*\[p = fS^i^) - ~S^{u*)\\p < a - ||5 - u*\\p. 
Such situation could happen only if Iju — u*\\p = 0, or u = u*. 

We can also receive an estimation ()4.9p from estimation (|4.10p as a limiting result for ^ — > oo. 

Now let us show that continuous dependence of operator S on t leads to continuous dependence of the 

=v . . 

solution of the problem on t. We will say that matrix integral operator S is continuous on t in the point 
to G [0,00). if for any sequence {tn} G [0;Oo) such as i„ — > to for n — > cxd. the following is true: 

(4.11) C(^) ^ 
for any u G 

From Theorem 2 it follows that for any t S [0,oo) equation (|4.4p has a unique solution, which is de- 
pending on t. Let us denote it as u^. We will say that solution of equation (|4.4p is continuously depending 
on t for t to, if for any sequence {t„} £ [0,oo); such as t„ to for n — >■ 00, the following is true: 



Theorem 3. Continuous dependence of solution on t. |21j 

Let us consider operator that satisfies condition (|4.7p for any t e [0,oo), with a independent on t 
-V 

and that operator Sf is continuous on t in a point to G [0,oo). Then for t ~ to the solution of equation 

(|4.4[) is continuously depending on t. 

Proof: 

Let us consider any t G [0,oo). We will construct the solution of equation (|4.4p as a limit of sequence 
{u„} : 



(4.12) 



— V 

Un+i = St (u„) (n = 0, 1,...; uq = 
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Let us rewrite inequality (|4.9p for n = 0: 

(4-13) - uqWp < ^ \\ui - uollp 

(1 - a) 

Since u^^ = Sf.^{Uf^), then because of (|4.12p and (|4.13p we have: 

(4.14) \\u; - u;x < (T^ll"! - "oll^' = (r^"^*'^"*"^ " ^^(^*o)llp 

Now with the help of (|4.1ip we obtain the required continuity of Ut for t = Iq. 
Let us consider now a norm L2 in space 

1/2 



Then we can rewrite estimation (|3.94p in the following form: 

Mx,t)\\h = I^In^Hx,t)\'dx^ < C. 
From the other side || • \\l2 is a weaker norm than \\ ■ \\p. 

Solution of the problem exists in norm || • \\p, which means that |iM(a;,t)||p < 00. Hence there exists 
solution in norm || • H^j, and hence condition p.94p 



(^Jj^,\u{x,t)\^dx^ < C 



is true for any t G [0,oo). 



Appendix A. 

The Fourier integral can be stated in the forms: 

-1 fOC /"OO /"OC 

U{ji,j2,l3)=F[u{xuX2,X3)]^j^^ I I I u{xi,X2,X3)e '^^'^^+^-^'+^'^''>dx,dX2dX3 



00 J —00 J —00 

00 nOO nOO 



-I pOO pOO pOO 

(A.l) u{xuX2,X3) = j:^-^ I I I (7(71,72,73)6 -*(^^^^+^^^^+^^^^)ci7id72d7; 



(27r) 

The Laplace integral is usually stated in the following form: 



'3 

00 J —00 J —00 



1 



(A.2) U^{rj)=L[u{t)]= I u{t)e-^*dt u{t) ^ — j U'^{7j)e "^drj c> cq 







c+zoo 



(A.3) L[u {t)]=rjU^{rj) -u{0) 
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The convolution theorem A.l. [32], 133] 

If integrals 

/•oo poo 
^0 "'0 

absolutely converge by Rerj > a a, then U^{r]) = Uf{vi) Uf{'q) is Laplace transform of 



(A.4) U{t)= I Ul{t-T)u2{T)dT 

Jo 

Useful Laplace integral: 



(A.5) L[e "^'l = / e *dt = ^ {RcT] > r]k) 

Jo [V - Vk) 
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